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Transformational optics allow for a markedly enhanced control of the electromagnetic wave trajectories within 
metamaterials with interesting applications ranging from perfect lenses to invisibility cloaks, carpets, concen- 
trators and rotators. Here, we present a review of curved anisotropic heterogeneous meta-surfaces designed 
using the tool of transformational plasmonics, in order to achieve a similar control for surface plasmon polari- 
tons in cylindrical and conical carpets, as well as cylindrical cloaks, concentrators and rotators of a non-convex 
cross-section. Finally, we provide an asymptotic form of the geometric potential for surface plasmon polaritons 
on such surfaces in the limit of small curvature. 
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1. Introduction to transformational optics 

It has been known since the advent of mathematical optics (P) that lenses with certain curva- 
ture and/or a spatially varying refractive index allow for focussed images. Indeed, the simplest 
example of a convergent lens is a curved piece of silica which bends light according to the Snell- 
Descartes laws of refraction at curved air/glass interface. However, such a lens is necessarily 
of a convex shape. On the contrary, a medium with a continuously changing refractive index 
can act as a flat convergent lens: Gradient-index optics is indeed the branch of optics covering 
optical effects produced by a gradual variation of the refractive index of a material with the 
famous paradigm of the Maxwell fisheye, first published by James Clerck Maxwell in 1854 under 
the pseudonym of Charles Anthony Swainson {2), but however still of current interest (fi,Q). A 
continuously refracting index can also create optical illusions known as mirages: in the desert, 
these are due to the variation of refractive index between the hot, less dense air at the surface 
of the road, and the denser cool air above it. The variation in temperature (and thus density) 
of the air actually causes a gradient in its refractive index, causing it to increase with height. 
This index gradient causes anomalous refraction of light rays (at a shallow angle to the road) 
from the sky, bending them into eye of the viewer, with their apparent location being the road's 
surface. Similar effects occur at the surface of the sea, and are known as Fata Morgana, after the 
name of the fairy in the Arthurian legend who lifted castles over lakes. Inverted optics include 
perfect lenses via negative refraction (5-:^), that can also be designed using transformational 
optics by mapping the image plane onto the source plane, which causes a negative (possibly spa- 
tially varying) refractive index (0, H)- Maybe the most stricking example of an optical illusion 
is the invisibility cloak independently proposed by Pendry, Schurig and Smith, and Leonhardt 
(fj^ . which can steer light around an object which is then invisible to an external observer 
(however, light emanating from this object is itself confined within the invisibility cloak!). The 
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underlying mechanism of such cloaks relies upon a mirage effect, since a source located inside 
the coating seems to radiate from a shifted location [Ph . 

This coating actually consists of a meta-material whose physical properties are deduced from 
a coordinate transformation in the Maxwell system, and this viewpoint allows for a refreshed 
analysis of geodesies of light in complex media, since the Maxwell equations retain their form 
under geometric changes and the only thing to do in the transformed coordinates is to replace the 
original medium (often isotropic homogeneous) by an anisotropic heterogeneous medium (73, U) 



a meta-material). Such geometric transforms are very useful in computational electromagnetism 



to model unbounded alS ) and twisted domains in the context of optical fibres. However, 
the physical interest in transformational optics has been fuelled by potential applications in 
the design of meta- materials, while the electrical engineering community has used geometric 
transforms as a mean to simplify numerical models involving large or thin domains for the past 
thirty years. 

Coming back to the invisibility cloak, the anisotropy and the heterogeneity of the parameters 
of its coat work as a deformation of the optical space around the object. The first experimental 
validation fj id ) of these theoretical considerations was given, a few months after the proposals 
of Pendry et al., and Leonhardt were released in the press. A copper cylinder was indeed made 
invisible to an incident plane wave at 8.5 GHz as predicted by the numerical simulations. This 
markedly enhance our capabilities to manipulate light, even in the extreme near field limit (f7^ ). 
However, such cloaks suffer from an inherent narrow bandwidth as their transformation optics 
design leads to singular tensors on the frontier of the invisibility region which is a consequence 
of tearing apart the metric when one makes a hole in optical space (obtained by blowing up a 
point ilS)). Physically, light has to curve its trajectory around the hole (or 'invisibility region'); 
Hence, to match the phase of a wave propagating in homogeneous space, it must travel faster. To 
remove the cloak's singularity, Xiang et al. proposed to consider the blowup of a segment instead 
of a point (f^ ). but this cloak only works for certain directions. On the other hand, Leonhardt 
and T yc c onsidered a stereographic projection of a virtual hyper-sphere in a four dimensional 
space (ia), which bears some resemblance with the construction of a Maxwell fisheye. 

As an alternative to non-singular cloaking Li and Pendry proposed a one-to-one geometric 
transform from a fiat to a curved ground: their invisibility carpet is by essence non singular 
and thus broadband. This proposal led to a rapid experimental progress in the construction of 



two-dimensional carpets approaching optical frequencies 



Another way to make cloaks broadband is to approximate their parameters using a homoge- 
nization approach, which leads to nearly ideal cloaking (j^ - [i^ ). as it does not rely upon locally 
resonant elements. In 2008, some of us 0) demonstrlTed broadband cloaking of transverse 
electric waves using a micro-structured metallic cloak in the micro-wave regime. This naturally 
prompts the question of whether, at optical frequencies, an object lying onto a metal film could 
be cloaked from a propagating Surface Plasmon Polaritons (SPPs). 



2. Prom transformational optics to plasmonics 

The less than usual physics of the transmission of light through holes small compared with the 
wavelength is has given rise to a vast amount of literature following the 1998 paper by Ebbesen et 
al. (fgil ) , but some heralding earlier work combining both theory and experiment is less well known 
(j^l ). Martin-Moreno and Garcia- Vidal further showed in 2004 that one can manipulate surface 
plasmon ad libitum via homogenization of structured surfaces (fg^j ). In the same vein, pioneering 
approaches to invisibility relying upon plasmonic metamaterials have already led to fascinating 
results (3^-32). These include plasmonic shells with a suitable out-of-phase polarizability in 



order to compensate the scattering from the knowledge of the electromagnetic parameters of the 
object to be hidden, and external cloaking, whereby a plasmonic resonance cancels the external 
field at the location of a set of electric dipoles. Recently, Baumeier et al. have demonstrated 
theoretically and experimentally that it is possible to reduce significantly the scattering of an 
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object by a SPP, when it is surrounded by two concentric rings of point scatterers [s^. 

In the present paper, we would hke to extend the control of electromagnetic fields using tools 
of geometric transforms to the area of surface plasmons. To do this, we adopt the general form 
for the time-harmonic electric and magnetic field, E and H in exp{jujt). 



H{x,y,z,t) = H{x,y,z)e 



j.t (1) 



where uj denotes the angular frequency and t is the time variable. We then consider the time- 
harmonic Maxwell's equations in the initial space coordinates (x,y,z): 



' \/ X E = -jujiH 

V X H = jujgE 
' V.igE) = ^"^^ 
V.ifiH) = 

where the tensors of permittivity and permeability e and describe the original (possibly het- 
erogeneous anisotropic) medium. 

In the transformed coordinate system (u,v,w), the structure of the Maxwell equations is pre- 
served: 



' V' X £; 

V xH 
' V'.{^E) 
V.{£H) 

but we note that the space derivatives are now taken with respect to the transformed coordinates 
(hence, the transformed gradient V), and obviously the permittivity and permeability tensors 
and /x' should have new expressions. 

An elegant way to identify the tensors is to consider the Jacobian matrix J associated with 
the change of co-ordinate system: 



ju}£E 







(3) 



dx dx dx 
ax = -TT-du + -TT-dv + TT—dw 
Qu Qv aw 

dy = T^du + T^dv + -^^dw 

QU QV QW 

UZ UZ uZ 

dz = -TT-du + T— + TT—dw 
ou ov aw 



dx\ 


/ du 


dy 


= jidv 


dz 


\ dw 



(4) 



It follows that dXi) the transformation rule for the tensors is: (e, fi) (e^, fj,') with: 



£!_{u,v,w) =J ^e{x,y,z)J -^det(J) 

fi'{u,v,w) = J^^fi{x,y,z)J~'^det(J) ^ ' 



where detf J) is the determinent of the Jacobian matrix. 
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When the original permittivity and permeabihty matrices are scalar, their transformed coun- 
terparts are given by: 

^ = eT'^ , ^ = fiT-^ , where T = -^-^ (6) 

We now wish to apply this general rule to the case of a transverse magnetic (p-polarized) SPP 
propagating in the positive x direction at the interface z = between metal {z < 0) and air 
(z>0): 

H2 = (0, Hy2,0) exp{i(A:a;2a; - ujt) - kz2z} ,z> , 

Hi = (0, Hyi,0) ex.p{i{kxix — ujt) + kziz} , z < , ^ ' 

where c is the speed of light in vacuum and £2 = 1 — u^+i-yui usual Drude form 

in the metal(z < 0), for which ojp is the plasma frequency (2175 THz) of the free electron 
gas and 7 is a characteristic collision frequency of about 4.35 THz (J^. Moreover, the upper 
medium is described by spatially varying tensors of permittivity and permeability /i' which 

for simplicity are assumed to be represented in a diagonal basis i.e. = d\a.g{exx2-,^yy2-,^zz2) 

and // = diag{Hxx2,fJ'yy2,fJ'Zz2)- 

The dispersion relation for the surface polariton at such an anisotropic interface takes the fol- 
lowing form (4% 



I £zz2£l(%y2gl 
C y El — exx2£zz2 

SPPs are bound to the interface, hence do not belong to the radiative spectrum (unlike leaky 
waves) . 



3. Design of plasmonic paradigms 

In the present section we show numerically that one can control the SPP field using transforma- 
tional plasmonics. A review of transformations leading to arbitrary carpets, concentrators and 
rotators is presented. 



3.1. Three-dimensional carpets for electromagnetic and plasmonic fields 

So far, such a mathematical derivation of the SPP dispersion relation seems fairly abstract. 
However, we now wish to analyse the interaction of this SPP with a specific anisotropic hetero- 
geneous structure, in the present case a three-dimensional invisibility carpet (j^), deduced from 
the following geometric transformation: 



y 

Z2 



(9) 



Z2 



Z + Z\. 



Importantly, we note that this material is not only heterogeneous anisotropic but also mag- 
netic, which seems a far technological reach. However, these constraints can be further relaxed 
using some quasi-conformal grids in the spirit of Li and Pendry's work for two-dimensional 
carpets ([H). 



November 10, 2010 1;17 Journal of Modern Optics mk8 



Journal of Modern Optics 5 

Here z' is a stretched coordinate. It is easily seen that this Unear geometric transform maps 
the surface zo{x,y) of the horizontal plan z{x,y) = onto the surface z{x,y) = zi{x,y), and it 
leaves the surface z{x,y) = Z2{x,y) unchanged. 

The surfaces zi and Z2 are assumed to be differentiable, and this ensures that the carpet won't 
display any singularity on its inner boundary. 

The symmetric tensors e[ and /x' are fully described by five non vanishing entries in a Cartesian 
basis: 



(a) Metal surface with a bump (b) Metal surface with a bump dressed with a 

transformed medium (invisibility carpet) 

Figure 1. Principle of a three-dimensional invisibility carpet for electromagnetic fields: (a) Light incident upon a curved 
surface undergoes different orders of diffraction (shown as black arrows); (b) Light incident upon a curved surface with an 
invisibility crapet (yellow box) undergoes the same diffraction as light incident upon a flat metal surface. 



Figure 2. Principle of a three-dimensional invisibility carpet for plasmonic fields: A surface plasmon polariton incident from 
the left (with trajectory shown as white arrows) on the metallic surface is smoothly bent in the metamaterial box (yellow 
box). It turns out that the transformed media in the yellow box is the same as in Figure [XT] 



/ 1 
1 

dz dz 1 



V dx' dy' 



a 







a 

dz dz 
y dx' dy' 



a 1 + 



dz 
dx' 



dz 

djc' 
dz 



+ 



idy') ) / 



(10) 



where a = {z2 — z\)l z^ and with J the Jacobian matrix of the transformation. Furthermore, the 
derivatives of z with respect to x' and y' are given by: 



dz 
dx' 



Z2 



Z2 dzi 



{z2 — z\)^ dx' ~^ ^ {z2 — zi)"^ dx' 



Zl — z' dZ2 



dz 
dy' 



Z2 



Z2 dzi Zl- z' dZ2 

7 + ^i7Tr-T\2 7^ (11) 



{z2 - z\f- dy' {z2 - z\f- dy' 
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The purpose of the next section is to show that such carpets work equally well for electromagnetic 
and plasmonic fields, as is intuitively the case when looking at the limit of ray optics, see Figures 
13.11 and 13.11 This is due to the fact that the transformed medium is valid for any field solution 
of the Maxwell's equations (which is obviously the case for SPPs). 



3.2. Gaussian shaped hump: 




Figure 3. Plane wave incident from the top (z axis): (a) Magnitude of the magnetic field; (b) Phase represention of the 
magnetic field. SPP launched from the left: (c) Slice of magnetic field along (x,y) i.e. top-view; (d) Slice of magnetic field 
along (x,z) i.e. side view. 



We now wish to apply the recipe for the design of 3D carpets to specific geometrical trans- 
formations in order to hide Gaussian and conical shaped bumps on a metallic surface. We first 
present the projection of the conical surface on the flat surface, as this case is more intuitive. 
The geometric transform used in such a design is as follows: 



/vrx\2 /7ry\ 

hn COS — — COS — — 



COS y— j with h 



2.10" 



/ = 1.25 10"^ Z2 = 7.5 10" 



(12) 



dzi 
dx 



= —2ho cos 



. sin 



/7rx\ vr 



(13) 



dzi [Try 



fTTX 

cos y— I • sm 



Try 



TT 
■J 



(14) 



We report in Figure [3^2] the finite element computations which exemplify the role played by the 
carpet in the control of the SPP wavefront. 



3.3. Conical shaped hump: A less than unusual carpet 

We now present the projection of the conical surface on the flat surface, as this case is less intu- 
itive. For this, we consider two surfaces zi and Z2 which describe the inner and outer boundaries 
of the carpet: 
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Figure 4. Plane wave incident from the top (z axis); (a) Magnitude of the magnetic field; (b) Phase represention of the 
magnetic field. SPP launched from the left: (c) Slice of magnetic field along (x,y) i.e. top-view; (d) Slice of magnetic field 
along (x,z) i.e. side view. 



zi = Aexp ( ] - b with ^ = 0.205 ; cr = 0.1 ; 6 = 0.005 



2 



Zi = h 



with h = OA; L= J2a^ln(^j^ 



(15) 



Exphcit expressions for the partial derivatives are 

dzi _ Ax f 1 + y2 \ dz-^^ ^ Ay / 1 + 

dz2 hx dz2 hy \ 

dx Ly/x^ + y2 ' dy L^J x'^ + y2 ' 

We report the finite element computations in Figure 13.31 We can see that control of electro- 
magnetic and plasmonic fields is achieved for both the smooth (Gaussian) and rough (conical 
shaped) bumps. 

We emphasize that these numerical results are non-intuitive as the scattering by a conical 
shaped surface is one of the hardest cases to handle. In order to see this, let us remind some 
classical result of singularity theory for a sector (the two-dimensional counterpart of a conical 
shaped domain). 

Let us consider the Laplace problem in a semi-infinite sector bounded by two edges Si and 

— An = / dans D , u = on Si and S2 , (17) 
which we can express in polar coordinates: 
d'^u 1 du 1 d'^u 

7^2 + -IT + = ^ Voni Q <r , Q<e <uj , (18) 

and let us assume that Dirichlet boundary conditions on the edges u(r, 0) = u{r,uj) = 0. 
Using separation of variables, we find that 



u[r, t) ) = 

n=l 



J2 (anr''^l"' + ftnr--""/'^) sin(n7r0/a;) . (19) 
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It is well-known (s^) that if / is square integrable, then u is of finite energy, that is 



giadu\^dx= I ( / ( (^) +(-^) ]de\rdr< +00 (20) 



/sector Jo \Jo \ dr r 89 

Therefore, we need to ensure that the following expressions are finite: 



-Y ( (an-r'"''/'^-^ - 6„-r-""/"-i)^^ rdr < +00 (21) 



^ V /"^ ^ f + 6nr-"-/-)') rdr < +00 . (22) 

2 „ Jo ^ ^ ^ 



We note that | En /o^^n(5^~"''^'^"^)^ ^'^^ is finite if and only if -Imr/u - 2 + 1 > -1. This 
condition is not met for < a; < 2tt. Hence, all coefficients bn vanish in (|19|) . 
The solution of the Dirichlet problem is thus given by: 



(r, e) = y] a„r""/'^ sin(— ) . (23) 

— ^ id 



n=0 



From this expression, we can see that the field u is singular near the tip of the sector i.e. when 
r is close to zero. Similar expresions hold for a conical shaped region in the three-dimensional 
case. However, when the sector /conical shaped region is filled with an heterogeneous anisotropic 
medium (for instance deduced from a geometric transform), the Laplace problem reads as: 

- V • (AVu) = / dans D , u = on Si and S2 , (24) 

where A is a matrix-valued function. The solution u is now very different from above, and its 
behaviour around the tip of the sector /conical shaped region therefore changes dramatically. 
This simple analysis exemplifies the role of the transformed medium in the regularization of the 
field scattered by a conical shaped bump surrounded by a carpet. 



3.4. Launching a SPP on a cylinder 

The goal of this section is to show that we can also launch a SPP on a cylindrical surface using 
the transformational optics approach in order to guide smoothly the SPP which propagates on 
a flat metal surface to the top of a cylindrical bump lying on the surface. 

More precisely, the transformed medium consists of two parts: On one hand, the quarter of a 
sphere (shown in pink color in figure [37S11 requires the following transform: 

zi = V0.04-x2 _ (0.22 -hy)2 ; Z2 = 0.5 (25) 

On the other hand, the cylindrical part (shown in yellow color in figure 13. 4p requires the 
following transform: 

zi = \/0.04 - 3;2 ; Z2 = 0.5 (26) 

We depict in Figure 13.41 the typical trajectory for a SPP launched on the tubular surface 
surrounded by the carpet consisting of these two parts. We report in Figure [3741 the associated 
finite element computations. 
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Figure 5. Principle of a hybrid tubular carpet for SPP. The transformed medium consists of two parts: The pink part 
amounts to projecting a semi-hemisphere on a flat surface, while the yellow part amounts to projecting one half of circular 
cylinder on a flat surface. The trajectory of the SPP is shown with white arrows for illustrative purpose. 




Figure 6. SPP incident from the left on the hybrid tubular bump of Figure [sTil which consists of one half of hemisphere 
and one half of circular cylinder lying on a metal surface. Tubular bump on its own (upper panel): (a) Slice of magnetic 
field along {x,z) i.e. side view; (b) Slice of magnetic field along (x,y) i.e. top view; Tubular bump with carpet (lower panel): 
(c) Slice of magnetic field along (x,z) i.e. side view; (d) Slice of magnetic field along (x,y) i.e. top view. 



3.5. Rotator 

In this section we consider an arbitrary shaped rotator described as follows: 

iii(6l) = 0.4 i? (1 + 0.2 sin(36')) ; i22(6') = (1 + 0.2 sin(36l) + 0.2 cos(46l)) ; ii = 0.4 . (27) 

One can see in Figure [3^ b) that the SPP is smoothly rotated within the rotator. Moreover, the 
concentrator is itself invisible for SPPs as the wavefronts are unperturbed outside this metama- 
terial. 

More precisely, we consider the geometric transform: 
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Figure 7. Principle of SPP scattering by a concentrator and a rotator: (a) The wavefront of an SPP incident from the top 
on a concentrator is squeezed within the meta-surface, while going unperturbed elsewhere; (b) The wavefront of an SPP 
incident from the left on a rotator is rotated within the meta-surface, while going unperturbed elsewhere. We note that in 
both systems, the wavefronts are parallel in the inner core, and are distorted inside the coating. 




m 



C22 



a.r 
1 + a^.r^ 



a.r 



V 



C22 





\ 



C22 / 



(29) 



C22 



de_ 



{R2{e) - Ri{e)) 



2 • 



dR2{e) dRi{9) 



86' 



86' 



{R2{6) - Ri{6)y 



R2{6) 



8Ri{6) 
86' 



Ri{6) 



dlM6y 
86' 



(30) 



3.6. Concentrator 

In this section we consider an arbitrary shaped concentrator described as fohows: 



Ri{6) = 0.4i?(l + 0.2sin(36')) ; R2{6) = 0.6i?(l + 0.2 sin(36l)) 
R = OA ■ Rs{e) = R{l + 0.2{sm{36)+cosii6))) 



(31) 



One can see in Figure [3T6l fa) that the SPP is smoothly squeezed within the concentrator. More- 
over, the concentrator is itself invisible for SPPs as the wavefronts are unperturbed outside this 
metamaterial. 

More precisely, we consider the following geometric transform: 



r' = a.r + (3 
a' — a 



with < 



a 



Ri{6) 



a 



R2{6) 
R3{6)-Ri{6) 

Rm-MG) 



(3 = 



(0 < r < Ri{6)) 



/5 = Me)^M^^ iRii6) < r < Rsi6)) 



R3i6) - R2{6) 



(32) 
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(a) Arbitrarily shaped concentrator (b) Arbitrarily shaped rotator 

Figure 8. Full-wave simulations for a concentrator and rotator: Phase representation of the y-polarized magnetic field for 
a SPP incident from the left (x-axis) with top view along (x,y) and side view along (x,z). (a) Concentrator; (b) Rotator. 



10 

Jrr' = I —a C22 
1 



/ (r - + 



,2 „,2 



m 



(r — /3).r 

C22-a 



r-P 
r 



r-l3 




/3 



/3 



R{ef (33) 



< r < Ri{e) : 



C22 



dr 

89' ^ ~Ri{e) 



09 



(34) 



Ri{e) < r < i?3(^) : 



+ - RimiRm - r) - ^^^(-R2(9) - -Ri(e))(iSi(9) - r)| (35) 



4. Geometric potential for SPPs on a weakly curved interface with an anisotropic 
medium 

In order to illustrate analytically the control of the evanescent component of the SPP, we consider 
the case of a curved surface which is invariant along the x-direction, see Figure |4l Following the 
paper by Delia Valle and Longhi (? ) where a detailed discussion of the role played by the 
geometric potential on the SPP envelope in the direction of propagation led to the evolution 
equation: 

OF d'^F 



in the case of two isotropic media. 
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Transformed medium (e2,M2) 




Metal (e^) 







Figure 9. Diagrammatic view of a metal- transformed medium interface and associated global and local coordinates systems. 



Delia Valle and Longhi used a two-scale asymptotic expansion in the electromagnetic field, 
assuming that the local curvature of the interface is small, i.e. | R ~ |S> 1 with R^^ = d(l)/da 
and a is a small positive parameter. Importanttly, the asymptotic expansions are assumed as 
follows: 

oo oo 

u = ^ a2i^(i) , V = ^ a2^+iv(2^+i) (37) 

i=0 j=0 

with u = {cBa, Eri, Ez) and v = {Ea,cBri,cBz) i.e. for p- and s-field components. Moreover, the 
corrective terms are functions of the slow variable cri = aa. 

The former authors obtain a hierarchy of auxilliary equations at successive powers of a. At 
the leading order (~ a^), they retrieve the usual effective index of the SPP: 



n 



eff = \/eie2/(ei + £2) (38) 



and proceed with the interface conditions for the field to obtain the asymptotic expansion of the 
decay length: 



7i = X;^?^^' = l'2 (39) 

i=0 



(7) 

where the corrective terms 7^ undergo some oscillations on the slow spatial scale cJi. 

Regarding the penetration depth of the SPP in the 2 media, Delia Valle and Longhi deduce 
the following expression: 

2 2e2,i - ej 2^181 
^^'^ = 2Rel-el ' ^^'^ 

from the solvability condition: 
at order a^. 

We performed a similar analysis for a transformed (anisotropic heterogeneous) medium above 
the curved surface. 

The transformed medium is assummed to be described by diagonal tensors of permittivity and 
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permeability 



§^{u,v,w) = diag(£(^,£^,£z) , //(u,v,w) = diag(/X(^, /x^, /x^) 



(42) 



where the geometric potential is given by: 




£l +£2 



1 



(43) 



In the local coordinate basis, we obtain the same evolution equation, with a new effective 



5. Concluding remarks and perspectives 

We have demonstrated the full control of surface plasmon polariton propagation by extend- 
ing the transformational optics tools to the area of plasmonics. We have shown that we can 
markedly reduce the scattering of SPP on a bumped surface thanks to invisibility carpets whose 
design is the same as for an out-of-plane electromagnetic field incident upon the metal surface. 
The versatility of the designs proposed in this paper illustrates the power of transformational 
plasmonics. 
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